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We have discovered from first-principles an unusual polar state in the low n Sr^+iTinOsn+i 
Ruddlesden-Popper (RP) layered perovskites in which ferroelectricity is nearly degenerate with 
antiferroelectricity, a relatively rare form of ferroic order. We show that epitaxial strain plays a 
key role in tuning the "perpendicular coherence length" of the ferroelectric mode, and does not 
induce ferroelectricity in these low dimensional RP materials as is well known to occur in SrTiOs. 
These systems present an opportunity to manipulate the coherence length of a ferroic distortion in 
a controlled way, without disorder or a free surface. 



Control over the emergence of (anti) ferroelectric and 
antiferrodistortive order remains a fundamental challenge 
for the atomic-scale rational design of new phenomena. 
Complex oxide heterostructures, layered thin- films, and 
other low-dimensional systems provide a novel platform 
to address this ongoing challenge. There have been two 
well-explored approaches to control ferroicity in these 
systems: epitaxial strain engineering - which has been 
used to induce ferroelectricity [1 , multiferroicity [2] and 
to create strongly coupled multiferroics [3] - and tailor- 
ing electrostatic boundary conditions. A significant chal- 
lenge in the field of ferroelectric thin-films has been un- 
derstanding the evolution of the spontaneous polarization 
(Pi in Fig. [T]) in ferroelectric/paraelectric heterostruc- 
tures [4] and ferroelectric nanoscale thin film capacitors 
as the dimension parallel to the polarization direction is 
reduced [SHS] • Today it is clear that this is driven almost 
entirely by the electrostatic depolarizing field arising at 
the interface. 

(Anti) Ferroelectric and antiferrodistortive distortions 
are cooperative phenomena involving the coherent mo- 
tion of atoms extending over many unit cells [10^ . There- 
fore, these ferroic states may also be greatly influenced by 
introducing a perturbation with a characteristic length 
scale below (or near) that of the coherence length of the 
ferroic distortion [11] [12]. Design strategies based on such 
an effect have lead to the emergence of many novel ferroic 
states, such as multiferroicity [13 and relaxor ferroelec- 
tricity. The ability to control ferroic order on the scale 
of the coherence length therefore represents an opportu- 
nity to understand and create novel ferroic phenomena. 
However, in bulk materials, the only known tuning mech- 
anisms are atomic disorder (which in many cases tends to 
introduce electronic defects that are detrimental to the 
desired properties) or free surfaces [M] . Clearly, alterna- 
tive pathways are desired. 

In this Letter we show how a complex oxide interface 
can exploit the coherence length of a ferroic instability 
in a controlled way and use it to design an unusual polar 
state in which ferroelectricity is nearly degenerate with 
antiferroelectricity, a relatively rare form of ferroic or- 
der. In contrast to the well studied "finite-size effects" 
problem in (anti)ferroelectrics [15 , the physics discussed 



here involves the emergence of local polar instabilities in 
a direction parallel to the interface. For such an in-plane 
direction, the polarization (P|| in Fig.jlj) never comes "in 
contact" with the interface, electrostatic boundary con- 
ditions correspond to a short circuit, perfect screening of 
the depolarization field is always achieved and the system 
is structurally infinite. 

As a model system we take the perovskite/rocksalt 
interface in the naturally occurring Ruddlesden-Popper 
(RP) layered perovskite, Sr^+iTinOsn+i- Dielectric 
studies and previous first-principles calculations have 
shown that low-n members of the series are paraelectric 
with low dielectric constants [161 This is surpris- 

ing given that the structure of the RP homologous series 
[TSl [TqI can be thought of as that of a SrTiOa perovskite 
with only an extra SrO rocksalt layer inserted every n 
perovskite unit cells along [001] as shown in the inset of 
Fig. [l] SrTiOa (the n = oo member of the RP series) 
is a well-known quantum paraelectric (QP) that displays 
a large dielectric constant (^^10,000 at low temperature) 
and can be driven ferroelectric with the application of a 
modest amount of biaxial strain ^ . 

The key discovery we make is that structural re- 
laxations occurring at the RP fault, that is, at the 
SrO/SrTiOa interface, break the coherence of the in- 
finitely long Ti-O-Ti chains parallel to the interface in 
different perovskite slabs. Even at a tensile strain value 
more than sufficient to drive SrTiOs ferroelectric (where 
the polarization lies in-plane, P||), here we find that the 
n = 1 RP is still very far from displaying a polar instabil- 
ity, even though the electrical boundary conditions and 
the length of the Ti-O-Ti chains in the relevant direction 
are identical to that of SrTiOs. An in-plane polar state 
must emerge with increasing n, the nature of which and 
how or why this happens is unknown. In the remainder 
of this Letter we explore these questions and show how, 
unlike in SrTiOs, epitaxial strain does not induce ferro- 
electricity in these materials. Instead, we elucidate the 
novel role of strain in tuning the perpendicular coherence 
length of the polar mode and use it to tune a system to 
a region of the phase diagram where ferroelectricity and 
antiferroelectricity compete with each other. 

We performed density-functional theory (DFT) calcu- 



2 



20000 




-5000 I ^ ^ ^ ^ ^ ^ ^ ^ ^ 1 

2 4 6 8 10 12 14 16 18 20 

n (# of perovskite blocks) 

FIG. 1: (Inset) The Ruddlesden- Popper, Sr^+iTinOsn+i, ho- 
mologous series showing two possible polarization directions: 
in-plane, Py, and out-of-plane, P±. (Main figure) Calcu- 
lated in-plane polar phonon frequency as a function of n for 
a fixed in-plane lattice constant of SrTiOs. (Blue dots are 
first-principles calculations, red curve is a fit to a single ex- 
ponential.) 



lations within the PBEsol approximation using PAW po- 
tentials, as implemented in VASP ^2]T[24] . The wave- 
functions were expanded in plane waves up to a cutoff of 
500 eV. Integrals over the Brillouin zone were approxi- 
mated by sums on a F-centered /c-point mesh consistent 
with an 8 X 8 X 8 mesh for the primitive perovskite unit 
cell. A low force threshold of 0.5 me v/A was used for 
all geometric relaxations in order to resolve the small en- 
ergy differences. Phonon frequencies and eigendisplace- 
ments were calculated using two methods: the direct 
method using symmetry adapted modes in VASP and 
Density Functional Perturbation Theory (DFPT) as im- 
plemented in the Quantum Espresso package. For the 
DFPT calculations, Vanderbilt Ultrasoft Pseudopoten- 
tials were used within Local Density Approximation. We 
ignore quantum fluctuations of nuclei and therefore bulk 
SrTiOa is predicted to have a ferroelectric ground state 
(See EPAPS). 

To begin to unravel the novel polar state that emerges 
in strained RP phases, we calculate the in-plane polar 
{Eu) phonon frequencies of the n = 1 to n = 5 members 
and the n = oo bulk. We fix the in-plane lattice con- 
stant of all RP structures to that of theoretical SrTiOs, 
a = 3.899 A. Note that the equilibrium in-plane lattice 
constant a increases monotonically with increasing n; the 
n = 1 member, Sr2Ti04, has a 0.5% smaller lattice con- 
stant than bulk perovskite SrTiOs. As Figure [l] shows, 
the soft mode frequency decreases monotonically with 
increasing n, even though the in-plane lattice constant 
was fixed to the same value for all n (and therefore the 
in-plane strain is actually decreasing with increasing n). 
Such a trend is surprising. The observed decay of the 



phonon frequency with increasing n can in fact be eas- 
ily modeled. It is exactly what one would find in a toy 
model calculation of the interplanar force constants of a 
finite thick, infinite slab as more layers are added to the 
slab (see EPAPS). It is, however, not clear why the RP 
phases, which are bulk materials, should display this kind 
of crossover from two-dimensional to three-dimensional 
behavior as n increases, given that the Ti-O-Ti chains 
are continuous and infinite parallel to the direction of 
the polar mode. Additionally, why does the polar mode 
of the n = 1 member have such a high frequency? We 
propose that in SrTiOs there is a coherence length per- 
pendicular to the direction of the polar mode and that 
coherence between different perovskite blocks is broken 
by the double rocksalt layers in the RP phases, effec- 
tively reducing the dimensionality of the system. There 
are therefore two questions that need to be answered: 
(1) does it make sense that a perpendicular coherence 
length exists in SrTiOs and if so can it be manipulated? 
and (2) can the double rocksalt layer really suppress the 
coherence between perovskite blocks? 

Critical thickness - Our basic premise is that the real- 
space coherence requirements of the lattice instabilities 
in the RP phases can be deduced from the phonon disper- 
sion curves of the bulk cubic perovskite. Figure [2^ shows 
that there are two main instabilities for bulk SrTiOs: an 
i?-point instability involving rotations of the oxygen octa- 
hedra, and the unstable polar mode at F. An alternative 
way to visualize the distribution of lattice instabilities in 
reciprocal space is by plotting the cj^ = isosurface in 
the first Brillouin zone [l2l|25], Fig.[2|b). The enclosed 
volumes on the zone boundary correspond to the i?-point 
octahedral rotation instability and the EE instability is 
associated with the volume in the zone center. We point 
out that this picture is similar to that obtained by Lasota 
et al, [26 who previously considered the coherence prop- 
erties of the i?-point rotation instability only. Fig. [2|a) 
shows that the branch stemming from the F instability 
becomes stable at wave vectors away from F in all the of 
high symmetry directions - towards X, M or R - and is 
therefore localized in reciprocal space to a finite volume 
around the zone center. 

From Fig. ^h) it is seen that the volume of the EE in- 
stability has a strongly anisotropic structure: it consists 
of three perpendicular disks. Each disk encloses wave 
vectors that have an instability involving ionic displace- 
ments in the direction perpendicular to the plane of the 
disk. The finite thickness of the disks, corresponding to a 
finite longitudinal coherence length, /|| ~ has been 

previously discussed in perovskite ferroelectrics such as 
BaTiOs [E] and KNbOs ES]. The finite radius of the un- 
stable disks in the reciprocal space found in SrTiOs sug- 
gests that the frequency of the unstable branch also de- 
pends on the perpendicular component of the wavevector, 
corresponding to a finite perpendicular coherence length, 

- l/q^ (See EPAPS). In contrast, in BaTiOs [T2l and 
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FIG. 2: Phonon dispersions of SrTiOs in its cubic (PmSm) 
phase at (a) experimental volume, and under (b) 0.6 % and (c) 
1.7 % isotropic tensile strain with respect to the experimental 
lattice constant. The corresponding cj^ = isosurfaces are 
shown next to each phonon dispersion curve. 



KNbOa [25 the isosurfaces consist of three almost per- 
fectly flat slabs that are infinitely extended, indicating 
that the frequency of the unstable branch is the same 
regardless of the component of the wavevector that is 
perpendicular to the direction of the ionic displacements. 
This implies that the critical thickness for an infinite slab 
of BaTiOs is zero |27 . 

Assuming the rocksalt layers break the perpendicular 
coherence between different perovskite blocks (we will 
shortly prove this is true), the non- vanishing critical 
thickness required for ferroelectricity in SrTiOa imposes 
a geometric condition, ncrit, for RPs to display an in- 
plane ferroelectric instability; n < ncrit members cannot 
have an energy lowering FE distortion simply because 
the number of perovskite blocks between the rocksalt lay- 
ers do not satisfy the coherence condition. However, as 
n increases, the structures should get closer to the FE 
transition. This is exactly what we have shown in Fig. [l] 
where the soft mode frequency u vanishes at ncrit ^ 7, 
agreeing with the lower bound from the reciprocal space 
picture of the instability of SrTiOs. 

Tuning the coherence length - The effects of strain [28] 
and pressure [29 on ferroelectricity in perovskites are 
well known; they alter the balance between the short- 
range forces (which favor the centrosymmetric state) and 
long-range Coulomb forces (which favor the ferroelectric 
state). But what are their effects on the coherence re- 
quirements? In Fig. [2] (b) and (c) we plot the phonon dis- 
persion curves and the corresponding cj^ = isosurfaces 



of cubic SrTiOa for isotropic tensile strains of 0.6% and 
1.7%, resepctively. As shown, negative pressure not only 
softens the F instability (the well-known strain-induced 
ferroelectricity [20]) but also increases the volume of un- 
stable q-vectors in reciprocal space. By 0.6% strain, the 
instability reaches the X point and at 1.7% strain the 
branch is unstable and relatively dispersionless in the en- 
tire F-X-M plane. In fact, the o;^ = isosurface at this 
isotropic strain value resembles that of BaTiOs [T2] . 

These results suggest that under increasing biaxial, in- 
plane tensile strain [30 the critical thickness for in-plane 
ferroelectricity in SrTiOs, and hence ncrit in the RPs, 
should decrease and eventually vanish. This is clearly 
seen for the RPs in Fig. [sj^a). Note that at a tensile 
strain value of 1.1% - a strain more than sufficient to 
drive SrTiOs ferroelectric - the perpendicular coherence 
length in SrTiOs is still non-zero. This results in 
the n = 1 RP remaining very far from displaying a polar 
instability. It isn't until a strain of ~1.4% that lj_ ^ 
and the n = 1 structure develops a polar instability. 

As further proof of a crossover from 2d to 3d ferroelec- 
tric behavior (and the direct strain control thereof), we 
parameterize from first-principles a finite thick, infinite 
perovskite slab force constant model for which we add 
additional layers of SrO and Ti02 planes. (For details, 
see EPAPS.) The softest force constant eigenvalue is plot- 
ted in Fig. [3)3 for different strain values as the number of 
SrTiOs layers, n, is increased. Notice that the in-plane 
polar force constant in this 2d model slab (Fig. ^jp) are 
evolving with both n and strain exactly like the in-plane 
polar instabilities in the structurally 3d RPs (Fig. |3^). 
This is clear evidence of the proposed coherence physics 
and the direct control of with stain. (A similar calcula- 
tion for unstrained BaTiOs, which has a ~ 0, shows 
that a single Ti02 layer in bulk BaTiOs is unstable.) 

The effects of rumpling- These results suggest that 
the rocksalt layer in Sr^+iTi^Osn+i breaks the coher- 
ence between perovskite blocks. Why and how? Note 
that rumpling of the SrO layers is permitted by sym- 
metry in paraelectric RPs and that the distance along 
the [001] direction between the Sr and O ions within the 
rocksalt layer can be as large as 0.20A [31 and quickly 
gets smaller in the layers further away, Fig.gb). In RPs 
with n > ncrit, the in-plane polar displacements of Ti 
along the infinite Ti-O-Ti chains get smaller the closer 
the chain is to the rocksalt layer. We propose that the 
rumpling at the double SrO layers breaks the coherence 
between perovskite slabs. 

To test this hypothesis, we artificially zero the rum- 
pling by moving the Sr and Ti atoms to exactly the same 
plane as the oxygens, and repeat the phonon calculations. 
We find that even the n = 1 Sr2Ti04 under zero strain 
has an in-plane polar instability and it is in fact almost 
equal in magnitude to that of SrTiOs, that is, the decay 
of the phonon frequency with n disappears |32J . 

What are the implications of our findings? We propose 
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FIG. 3: (a) Ruddlesden-Popper, (SrTi03)n(SrO): In-plane 
polar soft mode frequency squared versus layering number 
n. (b) Perovskite Slab Model, (SrTi03)n: Lowest in- 
plane polar interplanar force constant matrix eigenvalue ver- 
sus thickness of SrTiOs. Strain values are given with respect 
to the lattice constant of SrTiOs. The lines are fits to expo- 
nentials. 



that strained RPs, which lack coherence across the rock- 
salt layers, can not be ferroelectric, but rather display a 
novel polar state. In Fig. ^ we plot the energy versus 
mode amplitude in, e.g., Sr4Ti307 (n = 3) for both the 
ferroelectric mode and an antiferroelectric distortion in- 
volving polar distortions of neighboring perovskite slabs 
antiparallel to each other (Fig. [4^ shows a comparison of 
these modes). Notice that the energy surfaces for these 
distortions are degenerate up to the precision of this plot. 
In Fig. [4]i we plot the results of a similar calculation ex- 
cept that the rumpling was artificially zeroed; without 
rumpling, the degeneracy of the AFE and FE modes is 
lifted. These calculations (we found similar results for 
n = 1, 2 and 3 at two different values of strain) prove 
that the RP fault indeed breaks the coherence between 
perovskite slabs, resulting not only in a suppression of 
polar distortions for n < ncrit, but also for n > ncrit 
polar distortions in neigboring perovskite slabs do not 
interact with each other significantly even at higher than 
quadratic order. All of this suggests that there are an in- 
finite number of degenerate states involving uncorrelated 
atomic-scale polar regions. We surmise that the ground 
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FIG. 4: Sr4Ti30io ground state under 1.1% tensile strain, (a) 
Schematic of ferroelectric (FE) and antiferroelectric (AFE) 
distortions, (b) Rumpling distortion (defined as the distance 
along [001] between cations and oxygens on the same layer), 
(c) Energy gain due to FE and AFE distortions in structures 
with rumpling, (d) The same as (c) except the rumpling 
distortion has been artificially set to zero. 



state may be a form of relaxor ferroelectricity without 
disorder, the consequences of which remain unclear. 
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AUXILIARY MATERIAL 

Perpendicular coherence length.- It is possible to 
obtain information about the perpendicular coherence 
properties of a ferroelectric instability from the phonon 
dispersions of the high symmetry (paraelectric) struc- 
ture. Let us consider a cubic structure with f , y and 
z axes aligned with the (100) directions, and focus on 



the z component of polarization only. Figure |5] includes 
the sketch of distortion patterns corresponding to various 
wave vectors with qz = 0. The argument for BaTiOs [1] 
or KNbOa [2] is the following: As the unstable branch is 
dispersionless in the whole F - X - M plane, all of the dis- 
tortion patterns in Fig. [sjare equally unstable. (See Fig. 
|6]for the labels of high symmetry points.) This indicates 
that whether a mode on the polar branch is unstable or 
not does not depend on the variation of displacement pat- 
tern in the direction perpendicular to the polarization, 
and so a one unit cell thick chain of atoms is unstable 
by itself - even if all the atoms in the rest of the crys- 
tal are fixed in their high symmetry positions. (One can 
imagine taking a superposition of phonon modes with dif- 
ferent wave vectors to obtain a distortion pattern where 
only ions on such a thin chain are displaced.) Hence, the 
perpendicular coherence lengths for BaTiOs and KNbOs 
are zero. In other words, in a BaTiOs or KNbOs crys- 
tal in its high symmetry paraelectric phase, a coherent 
displacement of atoms that form a single unit cell thick 
chain can reduce energy. 

The situation is different for (unstrained) SrTiOs- Re- 
ferring to Fig. 2a of the main Letter, while there is a F 
point instability with local polar atomic displacements as 
shown in Fig. 5(a)[ the branch stemming from this in- 
stability is stable at the X (Fig. |5(c)D and M (Fig. |5(e)D 
points. Therefore the polar instability is localized around 
the F-point in k- space, stiffening as one transverses about 
halfway along F - X (Fig. 5(b)) and also along F - M 
(Fig. |5(d)[ ) points. This indicates that the frequency 
of the polar branch is not independent of the perpen- 
dicular component of the wave vector. A superposition 
of modes over this restricted volume of unstable phonon 
modes results in a finite thickness of the Ti-O-Ti chains. 
The minimum thickness of this unstable distortion is in- 
versely proportional to the radius of the disks in figure 2a 
of main text {q± in figure [t]), as this radius sets the num- 
ber of unstable phonon modes that can be superposed to 
obtain a localized distortion. (The shortest wavelength 
of an unstable mode with = is ~ l/<7±.) Hence for 
unstrained SrTiOs to display a polar instability, a criti- 
cal thickness of unit cells in a direction perpendicular to 
the polarization direction, must be achieved. 

Note, however, that this thickness scale obtained 
from the reciprocal space picture of the instability is a 
lower limit to the critical thickness, as the reciprocal 
space unstable disks are not perfectly flat, and that the 
relevant phonon branch is not perfectly dispersionless. 

Interplanar force constants.- In this section we pa- 
rameterize from first-principles a perovskite slab force 
constant model (finite thickness along out-of-plane di- 
rection [001], infinite along in-plane directions [100] and 
[010]) for which we add additional layers of SrO and Ti02 
planes along [001] and calculate the critical thickness for 
SrTiOs to display a polar instability along an in-plane 
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{a)qx =0, qy = 




(b)gcc = 0, qy = 7r/2a 




{c)qx = 0, qy = Tv/a 




(d)qx = 7r/2a, qy = 7r/2a 




{e)qx = Tv/a, qy = tt / a 

FIG. 5: Displacement patterns for the branch stemming from 
the FE mode for various wave vectors. Each arrow represents 
the local polarization direction (direction which the cations 
are displaced) of a single primitive unit cell. 



direction. 

The energy of a crystal can be expanded, in terms of 
atomic displacements with respect to the high symmetry 
positions of ions, up to quadratic order, as 

= ^0 + E E E c'^"'' {mui, {R') (1) 

a,a' ^^^/ 

Here, Au^{R) is the displacement of the k}^ atom in the 
unit cell designated by the Bravais lattice vector in the 




FIG. 6: First Brillouin zone for bulk SrTiOs in space group 
Pm3m. 




FIG. 7: Projection of one of the disks in Figure 2a of main 
text on the Qx - Qz plane. 



a^^ direction, and C^^! {R^ R') are the interatomic force 
constants (IPC). In general, R = niai + n2a2 + nsas 
where are the primitive lattice vectors and are in- 
tegers. For a simple tetragonal system ai -L a2 -\- as- In 
order to be specific, let as be along [001], which is the 
out-of-plane direction. We consider in-plane distortions 
only, Au^{R) = 0, The high symmetry structures we con- 
sider have 4-fold rotational symmetry so we can assume 
Au'^{R) = without loss of generality. From now on we 
drop the superscript, so Au}.{R) = Au^{R). 

We next impose in-plane translational invariance: The 
displacement of an ion in the unit cell R = niEi -\-n2a2 + 
nsas does not depend on ni or n2, but only on 713. That 
is: 



Au^{R = niai + n2a2 + nsas) = Au^{ns) (2) 



7 



The energy can now be expressed as 



E = Ec. 



77,3,77,3 



^{ns,n'^)Aui{ns)Aui,{n'^) (3) 



where N is the number of unit cells in a layer, and 
C'/^m:'(^3, ^3) is the interplanar force constant per atom 
(IPFC for short). C^^'{ns,n^^) is related to the force on 
an ion when a whole layer of ions are displaced, and it is 
intensive in the thermodynamic limit. 

The matrix of IPFCs can be calculated using the direct 
method. For this purpose we built a 1 x 1 x 6 supercell of 
SrTiOa (Fig. |8^) extended in the [001] direction. In order 
to get all the necessary IPFCs; 5 different displacements 
need to be considered: 1 for Sr (Fig. [sJd), 1 for Ti (Fig. 
[s]^) and 3 for O's (Fig. [sji). For each of these, we did a 
DFT calculation and obtained a column in the table of 
IPFCs, presented in Table |l] for unstrained SrTiOs. 

Next step is to build the IPFC matrices for a certain 
number of layers; and diagonalize them to obtain the low- 
est force constant eigenvalue. Size of the IPFC matrix of 
a slab is determined by the number of atoms in unit area 
(one perovskite unit cell). For instance, the IPFC matrix 
corresponding to slab that consists of a single Ti02 layer 
is 3 X 3: 



(1) 



2.617 
-0.190 
0.304 



-0.190 0.304 
12.285 -5.453 
-5.453 4.115 



(4) 



Here the first column corresponds to the displacement 
of Ti ion, and the other two columns correspond to the 
displacements of the two O ions on the same layer. (All 
the components are given in units of eV/A^.) A two-layer 
slab, which consists of a Ti02 layer and a neighboring 
SrO layer, has a 5 x 5 IPFC matrix 



(2) 



2.617 
-0.190 

0.304 
-1.491 

0.119 



-0.190 
12.285 
-5.453 
-0.582 
2.672 



0.304 
-5.453 
4.115 
0.146 
0.445 



-1.491 
-0.582 
0.146 
3.483 
0.371 



0.119 \ 
-2.672 
0.445 
0.371 
4.824 



\ 0.119 -2.672 0.445 0.371 4.824 / 

(5) 

where the last two columns are for displacements of Sr 
and the O on the same layer. The smallest eigenvalues 
of these two matrices are 1.36 eV/A^ and 1.13 eV/A^ 
respectively. These values give the first two data points 
for 0.0% strain curve in Fig. 3b of the main text. 

Both (7(1) and (7(2) involve force constants between 
ions only on the same or neighboring atomic layers, in 
other words, in the same perovskite layer. (One per- 
ovskite layer is two atomic layers.) So, they can be ob- 
tained solely from the '0*^ layer' row in Table [ij The 
IPFC for a 3 atomic layer (Ti02 - SrO - Ti02) slab, on 
the other hand, involves force constants between different 
Ti02 layers as well: 



(3) 





2.617 


-0.190 


0.304 


-1.491 


0.119 


-0.025 


0.067 


-0.009 \ 




-0.190 


12.285 


-5.453 


-0.582 


-2.672 


0.067 


-0.126 


0.006 




0.304 


-5.453 


4.115 


0.146 


0.445 


-0.009 


0.006 


-0.090 




-1.491 


-0.582 


0.146 


3.483 


0.371 


-1.491 


-0.582 


0.146 




0.119 


-2.672 


0.445 


0.371 


4.824 


0.119 


-2.672 


0.445 




-0.025 


0.067 


-0.009 


-1.491 


0.119 


2.617 


-0.190 


0.304 




0.067 


-0.126 


0.006 


-0.582 


-2.672 


-0.190 


12.285 


-5.453 


V 


-0.009 


0.006 


-0.090 


0.146 


0.445 


0.304 


-5.453 


4.115 / 



(6) 



The last three columns/rows of this matrix are for the atoms on the second Ti02 layer, the corresponding off diagonal 
elements are taken from the '1^* layer' row in Table [ij Similarly, the 4 layer slab has two more column/rows: 



(4) 



/ 


2.617 


-0.190 


0.304 


-1.491 


0.119 


-0.025 


0.067 


-0.009 


0.000 


0.022 \ 




-0.190 


12.285 


-5.453 


-0.582 


-2.672 


0.067 


-0.126 


0.006 


-0.001 


-0.045 




0.304 


-5.453 


4.115 


0.146 


0.445 


-0.009 


0.006 


-0.090 


0.002 


0.007 




-1.491 


-0.582 


0.146 


3.483 


0.371 


-1.491 


-0.582 


0.146 


0.036 


-0.037 




0.119 


-2.672 


0.445 


0.371 


4.824 


0.119 


-2.672 


0.445 


-0.037 


-0.427 




-0.025 


0.067 


-0.009 


-1.491 


0.119 


2.617 


-0.190 


0.304 


-1.491 


0.119 




0.067 


-0.126 


0.006 


-0.582 


-2.672 


-0.190 


12.285 


-5.453 


-0.582 


-2.672 




-0.009 


0.006 


-0.090 


0.146 


0.445 


0.304 


-5.453 


4.115 


0.146 


0.445 




0.000 


-0.001 


0.002 


0.036 


-0.037 


-1.491 


-0.582 


0.146 


3.483 


0.371 


V 


0.022 


-0.045 


0.007 


-0.037 


-0.427 


0.119 


-2.672 


0.445 


0.371 


4.824 / 



(7) 



Note that the diagonal components, which correspond to the self force constants, repeat, while the off diagonal 
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(a) Undistorted 



(b) Sr Displacement 



^ 





• • • • • 



(c) Ti Displacement 



(d) O Displacement 



FIG. 8: The 1x1x6 perovskite supercell and 3 of the 5 in- 
equivalent displacements used to calculate the IPFCs. Green 
spheres correspond to Sr ions, blue ones to Ti ions, and red 
ones to O ions. 



components involve IPFCs from lower rows of Table [l| 
This way, by using Table [Tj it is possible to build the 
IPFC for a slab of arbitrary number of layers and obtain 
the other data points in Fig. 3b of the main text. 

Finally, this whole procedure should be repeated to 
obtain the IPFCs and the eigenvalues corresponding to 
different strain states, and hence the other curves in Fig. 
3b of the main text. 

Notes: The force constants presented in Table |l] have 
been symmetrized. This process ensures that all the C(^) 
matrices are perfectly symmetric (the symmetrization 
changed the lowest eigenvalues no more than ~ 1%). 
The slab thickness in Fig. 3b, given in number of SrTiOa 
layers, corresponds to half of the number of atomic SrO 
or Ti02 layers in the slab. Half integer values correspond 







Ti 






Sr 


Osr 






Ti 


2.617 


-0.190 


0.304 


-1.491 


0.119 








-0.190 


12.285 


-5.453 


-0.582 


-2.672 




layer 


OT^± 


0.304 


-5.453 


4.115 


0.146 


0.445 






Sr 


-1.491 


-0.582 


0.146 


3.483 


0.371 






Osr 


0.119 


-2.672 


0.445 


0.371 


4.824 






Ti 


-0.025 


0.067 


-0.009 


0.000 


0.022 








0.067 


-0.126 


0.006 


-0.001 


-0.045 




layer 




-0.009 


0.006 


-0.090 


0.002 


0.007 






Sr 


0.000 


-0.001 


0.002 


0.036 


-0.037 






Osr 


0.022 


-0.045 


0.007 


-0.037 


-0.427 






Ti 


-0.005 


0.005 


-0.001 


0.001 


0.003 








0.005 


-0.005 


0.001 


0.001 


-0.003 


2nd 


layer 




-0.001 


0.001 


-0.001 


0.001 


0.000 






Sr 


0.001 


0.001 


0.001 


0.000 


-0.001 






Osr 


0.003 


-0.003 


0.000 


-0.001 


-0.019 






Ti 


-0.001 


0.001 


-0.001 












0.001 


0.000 


0.001 






3rd 


layer 




-0.001 


0.001 


0.000 










Sr 








0.001 


0.000 






Osr 








0.000 


-0.003 



TABLE I: Interplanar force constants for the unstrained 
structure, in units of eV/A^/ion. Columns denote the ion 
that is moved, rows denote the ion that the force acts on. 
Layers are listed in increasing distance from the displaced ion, 
and due to the finite size of the supercell, force constants for 
atomic layers as far as 3 unit cells are calculated. Subscripts 
for O's denote whether they are on a Ti02 or an SrO plane, 
and Oxiii denotes the force on the oxygen on the Ti-0 chain 
parallel to the displacement direction, Oxi^ denotes the other 
oxygen on the Ti02 layer. 
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FIG. 9: Energy gain, per formula unit, of the ferroelectric 
mode for n = 1 to n = 5 RP phases under different strain 
states. 



to an odd number of atomic layers that are terminated 
by Ti02 layers o both sides. Similar calculations for 
slabs terminated by SrO layers, not presented here, 
give qualitatively similar results even though the force 
constants are higher for small thicknesses. 

Transition Temperature estimate.- Fluctuations, 
both quantum [3] and thermal, are absence in DFT, so 
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the presence of a ferroelectric instability does not neces- 
sarily imply that the ground state is ferroelectric, even 
in the absence of other competing lattice instabilities. In 
order to obtain a crude estimate of transition tempera- 
tures, we present in figure [9] the energy gain of the FE 
state with respect to the PE state, as a function of n, 
calculated from first principles. 

Note that the energy gain for unstrained bulk SrTiOa 
is ~ 0.4 meV/f.u., and so we assume that quantum fiuc- 
tuations are more than enough to suppress an instability 
with this energy scale. This value is comparable to that 
of n = 1 under 1.7 % strain and n = 2 under 1.1 % strain. 
We therefore expect RPs with n < 2 and n < 3 not to 
distort in an experiment under 1.7 % strain and 1.1 % 
strain respectively. 

Another point of comparison is the energy gain 
calculated from first-principles, ~ 8 meV/f.u., for 
1.1% strained SrTiOa, which is ferroelectric at room 



temperature experimentally [4]. This is shown as the 
dotted line in figure [9| RPs with comparable energy 
gains, such as n > 4 under 1.1% strain or n > 2 under 
1.7 % strain are expected to undergo transitions on the 
order of room temperature. 
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